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Abstract 

A two-parametric family of integrable models (the SS model) that contains as particular cases several 
well known integrable quantum field theories is considered. After the quantum group restriction it 
describes a wide class of integrable perturbed conformal field theories. Exponential fields in the SS 
model are closely related to the primary fields in these perturbed theories. We use the bosonization 
approach to derive an integral representation for the form factors of the exponential fields in the SS 
model. The same representations for the sausage model and the cosine-cosine model are obtained as 
limiting cases. The results are tested at the special points, where the theory contains free particles. 

1. Introduction 

A complete set of form factors of all local and quasilocal operators, defined as matrix elements in the basis 
of asymptotic states, completely determines a model of quantum field theory and provides an important 
tool for studying its physical properties. In the case of an integrable massive model the form factors can 
be found exactly [1-3], as soon as the exact scattering matrix is known, as solutions of a set of functional 
(difference) equations named the form factor axioms. It is assumed (though not proven) that all solutions 
of these equations provide form factors of all local operators in the theory. 

There are several techniques for solving the form factor equations. In the case of a system of neutral 
particles of different masses, where the two-particle S matrices are scalar functions, any form factor factorizes 
into a product of given functions ('minimal form factors') and a trigonometric polynomial. The functional 
equations are reduced to an infinite chain of linear equations for coefficients of these polynomials [4]. For 
some particular local operators these equations can be solved for the whole set of form factors [5,6]. 

The situation is more complicated in the case of particles with isotopic degrees of freedom, where the S 
matrices are matrix functions while the form factors are multivector functions. In this case the form factor 
axioms have the form of equations in multidimensional spaces with matrix coefficients. The historically first 
approach to these systems was Smirnov's integral representation [3]. It gives form factors in terms of multiple 
integrals. The integrand can be represented as a particular transcendental function times a function of a 
given class, which depends on the local operator. In this way many integrable models have been studied. 
The advantage of this scheme is that it provides a general solution that corresponds to the most general 
local operator in the theory. The main problem of this approach is identification of sets of form factors with 
particular local operators in the model. It was only made for several most important operators like currents 
and the energy-momentum tensor. 

Another approach was proposed by Lukyanov [7] . It is based on the bosonization (free field representation) 
technique. Bosonization was first used in conformal field theory (CFT) as a way to construct an integral 
representation for correlation functions in the Virasoro minimal models [8] . Lukyanov showed that solutions 
to the form factor axioms can be found as traces of some objects, namely, vertex and screening operators, 
very similar to those of CFT. These operators can be expressed in terms of free bosonic operators. This 



provides integral representations for form factors different from those proposed by Smirnov. In the framework 
of this approach, the form factors of a wide class of local operators, namely the exponential fields, were found 
for the sine-Gordon model and the A series of the affine Toda models [9-11]. Note, that this approach is 
intimately related to the bosonization schemes proposed for lattice models of statistical mechanics [12]. 

It is worth to mention a new approach developed by Babujian and Karowski [13], based on what the 
authors call the off-shell Bethe ansatz. The advantage of this method is that is has to do with the same 
objects as the algebraic Bethe ansatz does. 

In this paper we construct an integral representation for the form factors of local fields in the two- 
parametric family of integrable field theories, which is also known as the SS model. This notation of 
the theory is related with the form of the scattering matrix for the fundamental particles. It is a tensor 
product of two sine-Gordon soliton S matrices with different coupling constants. The SS model has an 
explicit C/(l) (8" f/(l) symmetry, which can be extended up to the symmetry generated by two quantum affine 
algebras Uq-^ (sh) ® Uq^ (5^2)- This theory contains as the particular cases N — 2 supersymmetric sine-Gordon 
theory, 0(4) and 0(3) non-linear sigma models and several other interesting integrable theories. The SS 
model possesses a dual sigma model representation [15], which is useful for the short-distance renormalization 
group analysis of the theory. This analysis together with the conformal perturbation theory gives us a rather 
complete description of the short-distance properties of the model. The long-distance behavior can be derived 
from the form factor decomposition for correlation functions of the theory. 

The study of the form factors in the SS model was started by Smirnov [14], who calculated the matrix 
elements for a set of quantum group invariant operators as, for example, stress-energy tensor and U{1) cur- 
rents. An important class of operators in the SS model is formed by the exponential fields. In the restricted 
versions of the theory they correspond to primary fields in the related perturbed CFTs. Construction of an 
integral representation for the form factors of the exponential fields is the main problem that we consider in 
this paper. For this purpose we use a modification of the bosonization techniques proposed by Konno [20] 
for the analysis of integrable lattice models. 

The paper is organized as follows. In Sec. [5] we shortly give the necessary information about the model 
under consideration. In Sec.|31we describe the bosonization procedure for the model and construct a three- 
parameter family of form factors. In Sec^lwe identify the constructed form factors to those of the exponential 
operators. In Sees. El and we consider two degenerate limits of the model, namely, the sausage model and 
the cosine-cosine model, and give some checks of our construction. The computational details and some 
reference information are moved to appendices. 

2. The model 

The two-parametric family of integrable models under consideration (the SS model) possesses a Lagrangian 
formulation in terms of three scalar fields ipi, (p2, fs- The action has the form [15]: 

where the parameters ai, a2 and /3 satisfy the integrability condition 

al + al-p^ = 1. (2.2) 

It is convenient to introduce the notation 

Pi = 2al, P2^2al p3 = -2f3^, Pi + P2 + P3 ^ "2 (2.3) 

and 

as = -ip. (2.4) 

We have a two-parametric family of integrable models of quantum field theory. As particular cases this 
family covers a number of integrable families, like the sausage model [17] for p2 = 0, pi > 2, the cosine- 
cosine model [18] for pi + P2 = 2, pi,p2 > 0, the N = 2 supersymmetric sine-Gordon model [19] for 
P2 — 2 and pi > 0. Other models, like the parafermion sine-Gordon and an integrable perturbation of the 
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5{7(2)pj_2 X SU{2)p2-2/ SU{2)p-^-^-p2^4 coset model, can be obtained from this family by different restriction 
at appropriate values of pi, p2 [15]. 

There are three essentially different regimes in the theory: 

Pi,P2>0, P3 < (Regime I); (2.5a) 
Pi,P2,P3>0 (Regime II); (2.5b) 
pi,P2<0, P3>0 (Regime III). (2.5c) 

The regime I is unitary. Just this regime is the subject of this paper. The particle content and scattering 
theory in the nonunitary regimes II and III are rather complicated. The regime II is of particular interest, 
as it recovers the invariance under the substitutions ipi <-!■ ipj , pi pj . The spectrum and the S matrix of 
the model in the regime II, which can be extracted from the free field realization described below, is given 
in the Appendix IfI A detailed description of this regime will be given elsewhere. 

In the unitary regime the model possesses the U{1) x U{1) symmetry described by the conserved topo- 
logical charges^ 

Q± = i(Qi±Q2), Q.=Jdx'j°, ^^e'^-d,^,, i = l,2. (2.6) 

Classically the charges satisfy the conditions 

g±ez or Oi,Q2eZ, Qi + Q2e2Z. (2.7) 

In the quantum case these conditions are valid for the eigenvalues of the charges. Note that the operators Qi 
and Q2 are elements of a wider algebra, namely the above-mentioned Uq^{sl2) ® Uq^{sl2) with qi = — e'^^/P'. 

The spectrum of the model consists of the fundamental particles z^^i {e,e' = ± = ±1) and their bound 
states. The fundamental particles are characterized by eigenvalues of the topological charges: Qi\z^^i) = 
ejzee'), Q2\zee') = s'l^ee')- The mass of the fundamental quadruplet is proportional to the parameter fi of 
the Lagrangian. The exact relation between these quantities has a form [15]: 

„ r (21) r (2i) 

^ ^ M M 2 ; M 2 j ^ (2.8) 
The two-particle S matrix of the fundamental particles as a function of the rapidity difference 9 is given by 

Sp,pM = -SpAO)®SpM (2.9) 

with Sp{0) being the two-soliton S matrix of the sine-Gordon model with the coupling constant /3gQ — 
S^^^T [21]: 

siTiVi — 1 sm — 

S (0)++ - -e'^pW S (6)+- - -^e'^pi")— P— S (9)-+ - -e''^'^'^^) ^ 

p p 



r df 



, , , [°° dt sinh ^ sinh ZL(£_il£ sin 9t 



(2.10) 



sinh nt sinh 



The S matrix (|2.9|l was first considered in [14] and identified as the S matrix of the model H2.1|l in the 
region I in [15]. 

In the case when one of the parameters pi, p2 is less than one, the scattering matrix (|2.9f) possesses poles 
in the interior of the physical strip < Im 9 < n aX the points 

9 = \Ui^n, Ui,n ^ TT - npiu, n = 1,2,3,..., npi<l, if < pi < 1. (2-11) 



^In this paper we use the notation different from that used in [15]: the charges Q+ and Q— here are Qi and Q2 of the 
reference [15]. 
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These poles correspond to bound states. As in the unitary regime pi +P2 > 2, only one of these two sets of 
bound states can be present in the spectrum. The masses of the bound states are given by 

M,,„ = 2TOsin^!^, n = l,2,..., np, < I. (2.12) 

The bound states form quadruplets different from those of the fundamental particles. For example, for pi < 1 
each quadruplet consists of one particle with Q2 = 2, one particle with Q2 = —2 and two particles 
&g with Q2 — 0, while Qi = for the whole quadruplet. 



3. Bosonization and Integral Representation for Form Factors 

To compute form factors of local operators we would like to apply the bosonization procedure based on 
algebraic methods for solving the form factor equations. It was observed by Lukyanov [7] that solutions to 
the form factor axioms can be found in terms of representations of the corner Hamiltonian H and the vertex 
operators that satisfy the commutation relations of the extended Zamolodchikov-Faddeev (ZF) algebra. 
Namely, the form factors of a given local operator are identified with traces of products of the vertex 
operators with the operator exp(— 27r_ff) taken in a representation that depends on the particular form of 
the local operator.^ A crucial step in this approach is construction of a free field realization of the extended 
algebra. This makes it possible to compute the traces defining the form factors by application of the Wick 
theorem. 

Consider an integrable quantum field theory with the spectrum consisting of a set of particles {zj} with 
the charge conjugation matrix Cjj.^ Assume that scattering of two particles zj with the rapidity di and zj 
with the rapidity 62 into the particles zji and zj' is described by the S matrix S{di — 92)jj satisfying the 
unitarity relation, crossing symmetry condition, and the Yang-Baxter equation. Then the corresponding 
extended ZF algebra is generated by the vertex operators Zi{9), the corner Hamiltonian H, and the central 
elements $7/ with the defining relations 

Zj{9,)Zj{92) - 5(01 - 92yi'j'Zj,{92)Zi,{9i), (3.1a) 
/'./' 

Zj{9')Zj{9) = - . +0(1), 9'^9 + i7r (3.1b) 

(7 — (7 — ITT 

[H, Zi{9)] = i^Zi{9) - njZi{9). (3.1c) 

0.9 

If the scattering matrix for two particles zj and zj has a pole at the point 9 = luf j corresponding to a set 
K, of bound states of the same mass, then the vertex operators should satisfy some bootstrap conditions. 
Namely, the vertex operator for the particles zj^ with K d )C appear in the fusion of the operators Zj and 
Zj. This condition can be expressed algebraically as follows. Let u± be defined by the equations 

jc sin u+ 771/ 

M+ + U_ = Ujj, — = , 

smM_ rrij 

where mj, mj are masses of the particles z/, zj. Then we have the relation 

Zi{9' + iu+)Zj{9-iu_) = -^Y.^u^K^^)+0{l), 9' ^9 (3.1d) 

with some constants r|^. These constants can be found from compatibility with the conditions (|3.1a|) 
and Ij3.rbp and expressed in terms of the residues of the S matrix. Therefore, an important problem is to 
construct the basic vertex operators for the fundamental particles. Then the vertex operators for bound 
state particles can be derived by fusion of the basic ones. 

^This picture is deeply related [22] to the angular quantization procedure [23] where the operator exp(— 27ri?) is treated as 
a density matrix. 

^The charge conjugation matrix is related to the S matrix by the crossing symmetry condition: j,, C ji jn S(m — j = 

Y.j"C.j.r,s{eYj':i' . 
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Now we are in position to discuss application of the representation theory of the extended algebra 1)3. 1|) 
to form factors. For any representation 11 and any element X of the algebra we introduce the following 
notation 

((^))n = Trn(e-^-^) ■ ^^'^^ 

The main observation of the Ref. [7] was that the expression of the form {{Zij^{6n) ■ ■ ■ Zi^{9i)J)yy satisfies 
the form factor axioms, if we assume that it is a meromorphic function of the variables 0i, . . . ,9n- For 
9'i < . . . < O'j^-^, 6i < . . . < On the (m + n)-particlc form factor of some operator O is given by the relation 

= No{{Z''^ {9',+ '^)... Z''- (C + f )^/„(^n - f ) • ■ • Zi,(9^ ~ f )»n(o), (3-3) 



where 



Z\9) = (C-^y-^Zj{9). 



To obtain the form factors of a particular local operator O we have to choose appropriately the representation 
11(0) of the algebra It is important that in this formulation the form factor axioms in integrable 

quantum field theory are a direct consequence of the relations and the cyclic property of the trace. 

The normalization factor Nq in Eq. 1)3. 3|l cannot be established from the form factor axioms and should 
be calculated by other methods. The value of the central element J7/ on the representation Ifci gives the 
mutual locality index e^'^'^^ of the operator O and the local field V/(x) that creates the particle zj. It means 
that in the Euclidean space 

O(e2^'z,e-2'^'z)V/(0,0) = e2'^'^^C'(z,z)V/(0,0), z = .t^ + ix^, z^x^-\x^. (3.4) 

If the operator O has a fixed Lorentzian spin s the form factors possess an additional homogeneity property: 

iZi,(9N + -9)... Zi,(9r + 79))) - e^''((Z7„(0Ar) . . . Zi,(9r))). (3.5) 

For spinless operators (like the exponential operators, which will be considered below) the form factors are 
invariant with respect to overall shifts of rapidities and we have 

(vac|O(O)|ei/i,...,0jv/iv) = A^o((^/„(^Jv)...^7i(f?i)» for .s = 0. (3.6) 

The second important observation of the reference [7] is that the extended ZF algebras can be realized 
in terms of free bosons acting on a completion of a Fock space. Below we describe such a realization for the 
algebra H3.1|l of the model . 

In the case of the model 1)2.111 we have to introduce the vertex operators Z^e' iff) corresponding to the fun- 
damental particles Zgg' . The respective S matrix of the fundamental particles is Sp^p^ (9) given by Eqs. ()2.9)l . 
()2.1U)) and the charge conjugation matrix is 

For the values of the parameters, where the bound states appear, e. g. pi < 1, the matrix elements which 
contain the bound states corresponding to the poles ()2.11)l must be taken into account. The respective vertex 
operators can be expressed in terms of bilinear combinations of the operators for the fundamental particles: 

B'^"i9) ^ -irr.,\ Res Z++ie' + '-^)Z^+{9 (3.8a) 
Bl-"{9) - -i(i^s''"ri,„)"' Res{Z++{9' + %i)Z__(^? - %i) + Z+^{0' + %i)Z_+(0 - '-^)), (3.8b) 
bY{9) = -i(xi'"ri,„)-i Rcs{Z++{9' + %^)Z__(0 - %^) - Z+.{0' + %^)Z_+(0 - ii^)), (3.8c) 

u —6 

b]:'\9) = -irr,^ Res Z+^{9' + '^)Z^^{0 - %i) (3.8d) 
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with 

Tf^^i Res {^Sp,i9)-+SpAS)tX) = sin - • 5^, (mi,„)++5p, (mi,„)ii, 

/ ■ ,r , • «l „ \ 1/2 / • TT • «l „ \ 1/2 (3.9) 

/ sm — + sm — ^ \ , / sm sm — ^ \ ^ 



The nonzero elements of the charge conjugation matrix for any quadruplet of this type read 

Ch — = C_+ = Css = Caa = 1- 

In the case P2 < 1 the vertex operators B'^"{9), i3g'"(6'), B^"{0) are defined similarly. From the physical 
point of view both cases pi < 1 and P2 < 1 are equivalent, but the bosonization we develop below is not 
symmetric with respect to exchange of pi and p2- This is why we use both cases below while considering 
different degenerate limiting cases. The S matrices involving bound states can be found from Eqs. H3.1(l 
and 1)3. 8|l . We do not discuss them here in general because of their complicacy. Some important particular 
cases will be considered in Sees. El El 

Note that the vertex operators B^'" {ly = +, S, A, — ) given by Eq. I|3.8d|) is not the unique choice of vertex 
operators for the bound states. Linear transformations that do not affect the normalization condition 1)3. lb|) 
do not affect the A^-particle form factor contributions into the correlation functions. We shall use this fact 
to get rid of cumbersome coefficients in the vertex operators of the sausage model. 

Turn our attention to the operator contents of the theory (|2.1|l . A generic operator local with respect to 
the fields Lpi{x) is a linear combination of the monomials 

with arbitrary nonnegative integers kj. Finding form factors of all these operators is a difficult problem. 
Here we address a simpler problem of finding form factors for the exponential fields: 

To make the notations more symmetric it is convenient to introduce an imaginary parameter 

03 = -ih. (3.12) 

With this notation we have 

Oa,a,a,{x) = OjT/^, (x) = e-i<^i (-)+--'^- (-)+-3^3 (-) . (3.I3) 

The normalization factors for these operators were found in Rcf. [16]: 




, , ^ sr-^ smYi' aiait 2sinhPtn- sinh/'ii , o\ ,^ 

X exp / — — — > r-^ + ¥^ !--e~*> a- , 3.14) 

^' - 'sinht^ tanhafi sinh2(t/2) i J 

where 

r(z) 

P^^a^ai, Pi^P-2aia^, i = l,2,3; ^{z) ^ _ y 
i=i ^ ^' 

The products and sums are assumed to be over z = 1, 2, 3. 

Before describing the whole bosonization procedure, we give the resulting expressions for form factors. 
Consider a trace function {{Z^^^/^{9n) ■ ■ . .^e^e'^ (6'i)))jj(.0_^ „ „ )' which, according to Eqs. H3.3|l . 13.6|l . gives 
the form factors of the operator Oaia2a3{x) for N fundamental particles in the model H2.1|l . This function is 
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nonzero for even N with U{1) charges subject to the ncutrahty condition J^j = J2j ^'j = 0- F'o^' N — 2n 
we have 

A, A.. 1- 1 -J 1 J T)^ 



,Ari = ± s=l 



X 



where 6 = {9i, . . . , 02n), while the bold letters ^ = (^i, . . . , etc. are n-tuples. The values ki (and Ki which 
appear below) are defined by Eq. I|4.3|l . The numbers ji < . . . < jn {j[ < ■ ■ ■ < j'^) are defined so that 
Ej = — {e'j — — ) for j e {jsVs^i (j € {is}s=i) ^'^d — + Wj — +) otherwise. The functions in the r. h. s. 
read 

C = ciC2CiC2Cl (3.16) 

l<j<j'<2n 

= n A -/-iW2 n ^(^'^^^ - n ^^^^ - ^^■) > (^.is) 

n 

l<s<s' <n 

X n + (3.19) 

Here the signs ± are identified with ±1. The constants C3, Ci and functions G33((?), W(p]0), g\^^\9) can 
be found in Appendix IbI The constants Ci are given by Eq. H3.36() . The integration contours I?* = 'D\{A, B) 
in Eq. (|3.15() depend on summation term and are chosen as follows. The contour 'D\ goes from —00 to 
+00 below the poles of the integrand at the points Oj + '^{pi - 1 + 2Mpi) + 27riiV (Af, N = 0, 1, 2, . . .) 
and above the poles at 9j — ^(Pi ^ 1 + ^^Mpi) — 2tt\N . Note that the poles at 9j + ^{pi ~ + 2Mpi) are 
absent for j < js and the poles at 9j — ■^(Pi — 1 + 2Mpi) are absent for j > jg- The contour Vl goes 
above the points £,s' — Itt — 27riAi" and below the points ^s' + Itt + 27riAf, if Ag ~ Agi. A similar rule is valid 
for Pg. Besides, the contour goes below the poles at ^s' — ^{pi + P2) + 2mN and above the poles at 
^s' - y{pi + P2) - 2TTi{N + 1), if Bs = As' = +, and below the poles at ^s' + f{pi+P2) + ^mN and above 
the poles 6' + 'fiPi +P2) - 27ri{N + 1), if = A^' = -. 

Now we describe the construction of the form factors H3.15|) - (|3.19|) based on application of the bosoniza- 
tion procedure. Up to now, there is no general way to derive bosonization directly from the model and 
we need to guess it. So we introduce the construction in a mathematical manner, mostly without sub- 
stantiations. The definitions will be justified by the check of the relations (|3.1|l for the final boson vertex 
operators H3.35|l . As the construction is based on Konno's bosonization for the C/g,p(s^2) algebra [20], some 
motivations can be found in his paper. 

Consider the boson operators ai{t) (i Cz that depend on the real parameter t and satisfy the commu- 
tation relations 

sinh^ ^* 



am,a,{t')]^t . l^j,^, Sit + t')S,,. (3.20) 

smh vrr smh 



It is useful to introduce the fields 



MO\v)^ I ^a.(^)e'«*+-''l*l/^ (3.21a) 



— 00 



UO; ^a.(Oe'^*+-"l*l/^ (3.21b) 

J_oo 1^ smh^ 
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^<^\e; v)^2 sinh A,(±t)e^^<^*+-^*/\ (3.21c) 

We shall also use the notation 

xi^H^) = 0rH^;2-pO + '/'12(e;P.+i - 2) - -p,)- (3.21cl) 

In the expressions below we encounter the integrals 

dtfit) 

with f{t) having a pole at t = 0. We shall understand this integral as [24] 







with the contour Cq going from +oo + iO above the real axis, then around zero, and then below the real axis 
to +00 — iO. In fact, we shall use several simple formulas listed in the Appendix 1X1 
The Fock space is defined by the vacuum |0):^ 

a,{t)\0) = for t > 0. (3.22) 

as the space spanned by the vectors 

a.i^{-ti) . . .a.i^{~tn)\0), ii,...,i„>0, n = 0,1,2,.... 

This defines the evident normal ordering In what follows all traces are understood as traces over this 
Fock space. 

As building blocks for the generators of the algebra H3.1|l we introduce the operators 

Vi{9) = ■.exp{i(f>i+i{e;pi+i) + i(f)i+2{0; -Pi+2))-, (3.23a) 

ll^\e) = ■.eM~^M^■:P^) ± ^X^\d)):, (3.23b) 
These operators satisfy the following relations: 

V.{0')V,{9) - g,,{e - 9'y.V,{9)V,{e'y., (3.24a) 

v,{e')if\e) = wlf{e-e'):v,{e')if\0):, (3.24b) 

I^t^iO'me) - wlj\9 - e'):V,{e)I^^^\e'):, (3.24c) 

lr^\0')lf\e) = glf''\9 - e'):l[^\e')lf\e):. (3.24d) 

The functions gij{9) are defined as follows are understood modulo 3): 

dt sinh^ ^^ cosh ^ 



r°° rlf sinh 21 cosh ^ 

g^,[e) = G-\p,+,,0)G-\p.+2,e), G{p,e) = exp / ^ _^_^^e-^*, (3.25a) 

Jo t smh wt smh -f- 

g..,{B) = G-,\pk.e) {i^j, k + Gi{p,e) = exp / - l.pt ^~'''- (3-25b) 

Jo smh TTi smh 



dt sinh^ ^ 

The functions w^f '' (9) can be expressed in terms of the gamma-functions: 



4^^(^)=4r^W = l> (3-26a) 

4^Lm = wip^, m, wtlM = w{p,, m, (3.26b) 



* As our goal is calculation of traces, there is no need to introduce bra- vectors nor to discuss unitarity of the representations. 
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- = MP^, ^/m, (3.26c) 

where 



w{p,z\e) = rp = e(^^+i°g'^P)/P (3.26cl) 

with Ce bemg the Euler constant. Note, that all these functions have one series of poles at the points 
6 = —in + mpn or 9 = —in + inp{n + 1 /2) (n = 0, 1, 2, . . .) and one series of zeros at the points 6 = in — inpn 
or = Itt — inp(n + 1/2) . 

The functions g\f^\0) {A,B = ±) read 

sl'+'W =5(Pi,O,O|0), 5(r-)(^) =g(++)(0) ^ il5(p,,O,l|0), 5lr'(0) =.9(p„l,l|0), (3.27a) 

u — mpi^i/z 

(3.27b) 



with 



g(p,zi,Z2|g) = r| " (. (3.27c) 



r 



The functions g\'^^\{9) are rational and provide commutativity of the operators and /j^-* with i ^ j. It 
will be important later that the residues at the only poles of the products l[^\9)l[^\{9') satisfy the relation 

Res,,^,_^/i+H^?)4+^(^?') = -Res^,^^^^4-)(0)/i;)(0') (3.28) 
due to the important identity for the normal products 

■■l[''\0)l\tli(^~'-^)-- = ■.l\-\9)l[-l{9+'^):. (3.29) 

There is a general rule concerning positions of poles, which appear in the operator products. Let Uj{9) be 
normal ordered exponentials, like, for example, the operators Vi{9) or ij^^^ [9). Then the product Uk{9k)Uj{9j) 
can be reduced to the normal form 

Uk{9k)U,{9,) = gu.uASj - Ok) ■.Uu(9u)U,(9^):. 

The set of the poles of the function gUkUj i^j — Ok) in the variable 9j is bounded in the lower half plane, i. e. 
the imaginary parts of all poles are greater than some constant. On the contrary, the set of the poles in the 
variable 9k is bounded in the upper half plane. The same is true for any product Un{9n) ■ ■ ■ Ui{9i): the set 
of poles in the variable 9j arising due to any operator Uk{9k) with k > j (Uk is to the left of Uj) is bounded 
in the lower half plane, while the set of poles related to an operator Uk{9k) with k < j {Uk is to the right of 
Uj) is bounded in the upper half plane. 

As a result of the relations H3.24|l ~ (|3.27|l . the commutation relations of the operators introduced above 

are 

V.{0i)V,i92) = ~Sp^^,{9i - 92)XtSp^+M - 02)tXV,{92)V,i9i), (3.30a) 
V.i9i)V.±ii02) - ^'^^'^''l''l'h ±ii02)Vd9^), (3.30b) 



sinh 



92-6)i-i7r/2 



V.{9i)litli02) - . .,_,f+iW2 4+j(^2)V-.(gi), (3.30c) 



sinh ■ p 

cosh ^^-^^-W2 



Vd9i)ll^l{92) = ,,-,f+i./2 4^j(^2)V-(gi), (3.30d) 



cosh ■ 

p 
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J ^TTl 



sinh 



g-e-i7r/2 




cosh ^= — 

P2 



• e-if(p2-i) 

Fig. 1. Integration contours in the screening operators for different products that appear in Eqs. 13.3511 . 

sinh 

4^^(^i)4+i'(^2) = i\^l{e,)i\^\e,). (3.30f) 

Note that the commutation relations are completely independent of the values (±) of the superscripts. 

We can see from the relation H3.3()a|l that the operator V3 [6) is a good candidate for the operator Z++ (0) . 
However, it is known that the ZF algebras corresponding to nondiagonal S matrices cannot be represented in 
terms of pure exponentials of free fields. It is necessary to consider expressions containing integrals of normal 
exponents of free fields over the spectral variable 0. These integrals first appeared in CFT [8], where they 
screened the total charge in the auxiliary Coulomb gas and were called screening operators. In particular, in 
CFT they are necessary to provide the braiding relations for the vertex operators, which are CFT analogs 
of the ZF operators. 

In our case the analysis of the commutation relations for the vertex operators gives us the reasons to 
guess the following form for the screening operator: 

S,ik,n\9)^cJ -^(4+)(^)e«-i/i-)(0e-«)— ^^J^i— , (3.31) 

JCi ^"""i smh ^2 '— 

Pi 

with some normalization constants c^, which will be determined later. The contour Ci in this equation goes 
from —00 to +00 above the pole at the point 9 + i7r/2. As for the poles related to other operators, the 
contour goes below all poles arising due to the operators standing to the left of the screening operator Si 
and above the poles related to the operators standing to the right of Si. This is possible due to the remark 
concerning the position of the poles, which is formulated above after Eq. Ij3.26dp . 

In what follows the screening operators appear in the combinations Vi-i{9)Si{k, k\6) and Vi+i{6 + 
iTTpi/2)Si{k, k\9). In the first product the contour C; goes from —00 to +00 above the poles at the points 
9 + ^ — mpin {n = 0, 1, 2, . . .) and below those at the points 9 — ^ + mpin. In the second product it goes 
above the poles at 6* + ^ ~ nrpin and below those at 6* — ^ + ^T^Piin + 1). It is important to note that there 
is no real inflection of the contours between the poles at + i7r/2 and 9 — in/ 2 in the first product, because 
there is really only one of these two poles in each term: for /'•^■'(^) a pole may appear at 9 — i7r/2 and for 
at 9 + VK 12. In fact, the contour must be chosen separately for each term of the considered operator 
product (see Fig. Pl. 
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The meaning of the screening operator Si is to change the topological number Qi by —2.^ Each screening 
operator Si depend on its own pair of parameters ki and Hi. We shall see that, subject to some additional 
conditions, their values do not affect the commutation relations of the vertex operators. 

As the S matrix factorizes into a tensor product and the screening operators change the topological 
numbers, different screening operators must commute: 

[S,ih,K,\ei),Sj{k,,K,\92)]=0, iy^j. (3.32) 

This condition imposes some relations on the parameters ki, Ki. 

Since the screening currents l[^\9) commute according to Eq. (|3.30t|l . the only obstacle to commutativity 
of the screening operators is the poles of the functions g^it^\0) and g^[- \9). Due to these poles, we have 



[S'i(/Ci, Ki|6'i), S'2(/C2, K2\92) 



^ ^e''^+''"~''i^"''^(^""^^:/^^-'(^)4^-'(C - i^): 



7r^CiC2 



sinh«^^i^cosh«^^2^ 

Pl P2 



As, according to Eq. H3.29|l . the normal products in the right hand side coincide, the commutativity relation 
holds if 

To fulfill this equation we can take^ 

Kl + K2 = -^^2^2, 

and, by cyclic permutations, 

K2 + K3 = --fPsfca, + Hz = -^-piki. 

Solving these equations, we get 

ivr 

= -^{Piki +Pi+iki+i -pi+2ki+2)- (3.33) 

In what follows we assume Ki, K2, k.^ to be the functions of fci, k2, k^, given by H3.33|l . 

Now we are ready to express the vertex operators and the corner Hamiltonian in terms of the boson 
operators ai{t). To do this we introduce an auxiliary algebra generated by two elements lo and p with the 
relations 

up' — — 1, uop — —pijj, Trp = Troj = 0. (3.34) 
Then the generators of the algebra (|3.1|l for the model 1)2.1(1 can be represented in the form 

/ dtY, -^-^a,{-t)a,{t), (3.35a) 

Jo ~[ smh — 

Z++(fci,A;2,A:3|0) =c^V^3(^)e("^+'■^''/^ (3.35b) 

Z„+(fci,A;2,A:3|0) =c^pF3(0)5i(fci,«i|0)e(^i+'^^)'^/2, (3.35c) 

Z+_(fci,fe,fc3|0) = -,jpVmS2ik2,K2\9-'^Y^''-^''-^'l\ (3.35d) 

Z__(fci, fe, A:3|e) = -ujVz{,9)Sx{kr,Kr\9)S2{k2,K2\9 - il^)e(''^i+'=^)«/2. (3.35e) 

These operators satisfy the commutation relations ((3.1a(l and ((3.1c(l . In addition, if the normalization 
constants are given by 

g2(CE+log7rp,)/p. r(l + l/pi 



^3/2 r(-i/K) 



G(p„-i7r), (3.36) 



^The charge Q3 will appear in Sec.|l] when we will consider other regions of the parameters of the model 12.11 . 
® Other solutions do not lead to physically different results. 
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they satisfy the normahzation condition l)3.1b|l . The proofs of the commutation relations and normahzation 
condition is rather standard and are presented in Appendix lUl 

The expressions for the operators H3.35|l define a free field realization of the algebra H3.1|) for our model. 
This realization reduces the problem of construction of the integral representation for form factors to appli- 
cation of the Wick theorem. For arbitrary free field exponents (e. g. V{0), I^^\0)) of the form 



Uj{e) = :e^^(^': = e*^ ^^^e*^ ^f{e) ^ V / dt A){±t)a,{±t)e 

Jo 



±iet 



(3.37) 



the Wick theorem gives 



with 



{{UM . . . UM)) = n n ^u,u, (Ok - e,) (3.38) 

j=l j<k 



\ogCu,^{{cf>jm+{0))), (3.39a) 

logGu.uM^Ujmkm- (3.39b) 

The traces (I3.39f) can be easily calculated. The values of Cuj , GujUk ("") for the operators, introduced in the 
paper, are listed in the Appendix The expressions H3.15|l - (|3.19|l are obtained from these formulas after 
the substitution T i7rpi/2). These shifts of the integration variables allow one to deform 

two first contours depicted in Fig. ^ to the single straight contour like the third one in Fig. ^ 

We have a three-parametric family of form factors labeled by fci, fc2, k^. Now we have to identify them 
with a three-parametric family of local operators in the Lagrangian formulation. It will be done in the next 
section. 

4. Identification of form factors 



Ii 
h 
h 



Due to the formal symmetry of the action (|2.1|) with respect to the substitutions (ipi^ai) ^ (ipj^aj), the 
model is unitary in three regions: 

Pi < 0, P2,P3 > 0; 

P2<0, pi,p3>0; (4.1) 

P3 < 0, Pl,P2 > 0. 

What we considered above, while considering the regime I, was the region I3. 

In the region 1^ the operators corresponding to the fundamental particle is defined in terms of the operator 
Vi. Namely, we can introduce the vertex operators (recall that i S Z3) 

Zl+{h,k2,h\e) = t^F,(0)e('='+i+^'+^)«/2, (4.2a) 

^l+(fci, fc2, fcal^) = ^pl^.(0)5.+i(fc,+i, K,+i|0)e('='+i+'='+^)«/2, (4.2b) 

Z;_(fci,fc2,fc3|0) = -u;pV,{0)S,+2{h+2,n,+2\e - '^^)e^''^+'+'''+''>'/^ (4.2c) 

Zl_{ki,k2,ks\e) = -u;^,(0)5,+i(fc,+i,K,+i|0)5,+2(A:,+2,K.+2|e- ^^)e('='+^+"'+^)'/'. (4.2d) 

Evidently, in the previous section we constructed the generators Z^^,{ki,k2,k3\6). The vertex operators 
Zl^,{6) for a given i satisfy the relations H3.1|l with the 5 matrix —Sp^^^{9) (g) Sp^^^{9). The topological 
charges in these regions are given by the same formula H2.6|) for Qi, but without the limitation to i = 1, 2. 

We have a three-parametric family of form factors. We expect that they are in a one-to-one corre- 
spondence with the exponential operators H3.13|l . Our main conjecture is that the correspondence between 
oi, 0-2, ^3 and ki,k2, k^ analytically depends on the parameters pi, p2, Ps for any values of these parameters 
including all three regions Ii, I2, I3. This analyticity conjecture becomes more natural, if we consider the 
intermediate region of parameters, corresponding to the regime II, where all three families of vertex operators 
coexist (see Appendix IF|) . 
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Now we want to exploit the mutual locality indices, introduced in Eq. H3.4|l . For this purpose we need 
operators that create the fundamental particles z^^, . Such operators are given by 

v:,(.) = oL,(.)exp(i£|-iM + i£|^ 

with <fi{x) being the dual field to (pi{x), d^(pi{x) = e^'^ di,Lpi{x) and 0\^,{x) being an operator mutually local 
with the operators ^Pi{x). Indeed, consider the commutation relations [Q^, V*^/] in the Euclidean plane. We 
have 

[g, , VI,, (x)] = - ^ ^ dy^ d,^, (y) V:,, (x) . 

The last integral can be taken over a small circle around the point x, where we can neglect the interaction 
term in Eq. (|2.1|) and consider the model as that of three free bosons. It can be easily taken using the 
complex variables z = x^ + ix^, z — x^ — \x'^ with the result 

It means that this operator generally creates from the vacuum all states with topological charges Qi+i — s, 
Qi+2 = e' including the one-particle ones. 

The mutual locality index of the operator OaiaaogCa;) with the operator Vl,i{x) is equal to 

g27riOJ^, _ giTr(6ai+i/ai+i+e'ai+2/ai+2) 

On the other hand, using Eqs. Ij^.tcl) . it is possible to calculate the mutual locality index of the operator 
V*g/ (x) and the operator that corresponds to the form factors defined by the generators (|4.2I) with given ki , 
k2, fcs: 

g27rin^^, _ gi7r(efci+i+e'fci + 2) 

Hence, we conjecture that the form factors of the operator 0010203(2^) i^i the region I3 are given by the 
traces of the operators (|3.35|) with 

, a-t in f a, a,+i a,+2 \ 

ki ^ — , = -— Pi+2 ■ (4.3) 

tti A \ ai ai+i ai+2 J 



5. Sausage Model 

In the limit a2 ^ 0, the action 1)2. l|l splits into two parts 

S[(pi,ip2,ip3] = Ssm[vi,V3]+ / dx—^ — , (5.1) 

where Sqm is the action of the sausage model in the dual representation: 

Ssm[^,x]= I d'x (^ (^M^)' + (^mX)^ + ^cosa(^cosh/?x) , ~ ^ 1 (5.2) 

with a = ai. The mass term for the field 1^32(2;) in Eq. (|5.1|l appears in the limit P2 ^ due to the quantum 
corrections (see Ref. [15] for details). 

On the level of the factorized scattering theory the limit P2 ^ must be performed as follows. According 
to Eq. (|2.8(l . H2.12|l the mass of the fundamental particles m tends to infinity, while the masses of bound 
states remain finite. Hence, the fundamental particles decouple. The mass of the nth bound state is 2/Ltn, 
which means that nth bound state is simply n the first bound state particles with the same rapidity. It 
means that the model in this limit is completely described by the first bound state particles. 
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Consider the first bound state of the particles 2^^+ and z^^-- This bound state forms a quadruplet. Its 
scattering matrix is given by the diagram 




-P2) 



The values of the rapidities are written near the arrows. For P2 = the well known fusion phenomenon takes 
place. The quadruplet V'^ (8) splits into a triplet and a singlet V-^. The S matrix acts on the spaces 
^V^, ^V^, and (8) trivially and on the space as a. nontrivial S matrix. It means that 

the singlet ZQ{d) describes a free boson particle 1^2(2;) and can be decoupled. The remaining triplet Zi{9) 
(where the subscript / = +, 0, — labels the particles with the charge Qi = 2, 0, —2 respectively) corresponds 
to particles with the mass M related to the parameter /i in the action as 

M = 2fi. (5.3) 

Scattering of these particles is described by the S matrix of the sausage model SM\ [17] with 



11 , ^ 

X^- = —. 5.4 

p 2a^ 



The S matrix of the sausage model has the form [17,25] 
sinh^^ 

s{e)++ ^ 



+ sinh^' 

p 

p p 

sin ^ sin ^ sinh ^ sinh 

S{e)l+ = ^ V^'5(0)++ 8(9)+- = ^ ^^3(9)++, 

^ sinh ^ sinh ^ ^ sinh ^ sinh ^ 

p p p p 

i sin — sinh - 

= sinh^'sinh^ ^^^^++' = + 

p p 



(5.5) 



-I2 _ Q(0\Ill2 _ Q(0\l'll2 

and the charge conjugation matrix is 



C/,7 = d.ij. (5.6) 

To obtain the generators of the algebra (|3.1I) for the sausage model we have to consider the first bound 
state operators B^'^{9) [v — +,S,A, — ) for P2 < 1 and then take the limit p2 — ^ 0. The operators Bl "{9) 
can be obtained from the operators Bl'"{9) defined in Eq. H3.8|l by the substitutions pi ^ p2, Mi,n ~* W2,n 

and Zee'(^) ^ ^s's(^). 

In the limit P2 we easily obtain 



, ifg'^ - 2p2 ^/^./— sin— , A"^'^ ^ ./2cos— , (5.7) 
pisin^ \ Pi V Pi 



We see that the coefficients at Bj^ (9), Bj^ (9), and i?J (9) in equations analogous to Eq. (|3.8|l remain finite. 
These three operators form the triplet. The coefficients at the remaining operator Bg'^{9) is singular. This 
is just the singlet operator. 
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We note that the contributions of the oscillators ai{t) and 03 (i) to the vertex operators i?g'^ corresponding 
to the singlet particle vanish, and this operator can be written as 

/oo 
dta{t)e'^* +iV2^a2 (5.8) 
-00 

with the parameter 02 defined in Eq. (|3.13() and the operators a{t) related with the operators 02 (i) as 



a(0 ^ -|^^y "^'"^ ' a2it), [dit), dit' )] = 2 sinh TTt 6 (t + t'). (5.9) 

The corresponding traces consisting of the operators Zq{6) are given by 
((Z^(79„, + ■m/2) . . . Z'ai^i + i7r/2)Z^(^i - m/2) . . . 4(0„ - i7r/2))) 



= (1^2^02)"+" + ((5-functional terms). (5.10) 
They reproduce the form factors of a free field. Indeed, for a free boson field with the action 

'5M= / ^ ((9,^)^-MV^) (5.11) 



we have 



iniii(m,n) k 



k=0 l<ii<...<ife<™ 
i<ji<-<jk<" 



= (iV27r a)™+" + (J-functional terms). (5.12) 

It is easy to see that the (5-functional terms (which correspond to nonconnected diagrams) in Eq. (|5.1Q|I are 
the same as in Eq. H5.12|l . 

This provides a check for our construction. The consideration in Sec. 0] does not distinguish between 
primary (exponential) fields H3.13|l and zero spin linear combinations of their descendants (|3.10(l . The result 
(|5.10|l . as well as other free field cases described below, provide an evidence in favor of the identification of 
the form factors just with the primary operators. 

In the triplet operators B^-^{0) {v = +, A, — ) the contribution of the oscillators a2{t) vanishes, and we 
obtain the generators of the algebra H3.1|l for the sausage model in terms of two families of boson operators 
a{t) — ai{t) and b{t) = a3(<). 

Namely, let 



sinh 



2 TTt 
2 



ait),a{t')]^t———^—^S{t + t'), (5.13a) 



sinh TTt sinh 
sinh^ ^ 



2 



sinh TTt sinh 



The fields associated to the operators a{t) according to the rules H3.21|l will be denoted as (j), while those 
associated to the operators b{t) will be denoted as ip. Now we introduce some building blocks for the explicit 
representation of the triplet vertex operators. Let 

V^^\9) = :exp {i^{9;p) T i{4>^^^H0; 2 - p) ~ V^(±)(0; -p))):, (5.14a) 
I^^^\e) = ■.cxp{^i4>{d;p)±i{(P'^^\9;2^p)-^/^\9;-p))):. (5.14b) 

The operator products for these operators are given by 

V^^^e')V'^^H9)=g'^'^^\9-e'):V'^^\e')V^^H9y., (5.15a) 
^ (g'^^^'^e - 0')r^ ■.l'^^\0')v'-''\9y., (5.15b) 
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= .g(^^)(0 - e') ■.i^^\e')i^^\e):, (s.isd) 

with 

Then we define the 'bare' vertex and screening operators as 

V{k, k\B) = {v''+\e)e"' - iy(-)(6i)e'=) e'=^ (5.16) 

S±{k,n\e) = £ ^ (z'+HOe'' - i/^-HOe-'^) "^^i^- (5-17) 

The contour C+ (C_) goes below (above) the pole at the point 6 — iir (6* + Itt). 

The screening operators appear in the products 5+(fc, k\9)V''^^9) and V'^^^ {9)S^{k^ k\9). The products 
V^^\9)I^^^\S^) possess a simple pole at the point £, = 9. Taking this fact into account we should carefully 
define the prescription for the integration contours in the screening operators S±{k^ hi\9). According to the 
general rule, since the screening operator 5*+ {S-) is placed to the left (right) of the operator V'^^\9)^ the 
contour C+ (C_) goes above (below) the point 9. Besides, both of them go below the poles at the points 
— ivr + 'mpn {n = 0, 1, 2, . . .) and above the poles at + Itt — 'mpn. Similarly to the case of the operators 
5i, there arc no inflections in the contours, because the poles f = + Itt, — Itt never appear in the same 
product simultaneously. 

We can find from Eqs. (|5.15b|) . H5.15c() that 

V^^H9)l'^\i) . = } (5.18) 

smh ^= smh ^= ' — 

V p 

Using this relation we obtain 

S+{k, K\9)V{k, k\9) = V{k, n\9)S-{k, n\9), (5.19) 

The pole at = 6* is cancelled due to the identity g'^~^^^{9) — g^^^^ff). But for the product of the operator 
V{k,K\9) with two screening operators the pole contributions do not cancel. Namely, 

S+{k,K\9)(^S+{k,K\9)V{k,K\9) - V(k, K\e)S-{k, k\9)'^ = m^~''^''~'''U^+\9 ~ ni), 
(s+{k,n\9)V{k,n\9)-V{k,n\9)S-{k,K\9)^S-{k,K\9) = -i5e~'=-'=(''+^')/(-)(6' + ttI), (5-20) 

These relations take place due to pinching the contours of two screening operators between poles. The 
situation is analyzed in detail in Appendix IdI 

Now we are in position to introduce the generators of the algebra for the sausage model. They read 

Bl^{9)=c-^^Z+{k,K\9), Bl\e) = -Zo{k,K\9), B^_:\9) = ciZ^{k, k\9), (5.21) 
where the constant ci is defined in Eq. (|3.36() and the renormalized generators are 

/•°° /sinhTTtsinh^ , , , sinhTTtsinh^^^^^ , . , \ 

Z+(fc, ti\e) = cpV{k, Hi\9), (5.22b) 
Zo(fc, k\9) = -Cyj2cos^V{k, K\9)S-{k, k\9), (5.22c) 
Z_(fc, k\9) = ~-cpS+{k, K\9)V{k, K\9)S_{k, k\9). (5.22d) 
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Fig. 2. Integration contours in the terms l|5.2()|l consisting the operator Zo{9). The pictures are labeled 
by values of the pair (AB). On the left picture the solid line corresponds to the screening operator S-, 
while the dashed one to S+. 



As we have mentioned above, the definition of the bound state vertex operators is not unique. The operators 
Zj{k,K\9) (/ = +,0,—) can be considered as vertex operators for the sausage model, because the trans- 
formation H5.21|l respects the normalization condition (|3.1b|l . The guidelines for checking the ZF algebra 
commutation relations in this case are given in Appendix IdI 

The constants fc, k are related to the parameters of the exponential operators. For the operator 

0„fc(a;) =e'°'^(^)+''^(^) (5.23) 

we have 

k=-, n = ^-ip- + {p-2)-) . 5.24 
a 4 \ a p J 

The normalization constant c in Eq. H5.22II is given by 

g2(CB+log77p)/p rl/2(l/p) 



(7rp)3/2 r3/2(l-l/p)' 



(5.25) 



Now we make the integration contours for the operators Zq and Z_ more obvious. The operator Zq 
contains the terms of the form 

2m ^ ' ^^^sinhi^^ ^ ' 

p 

The integrations contours for these terms labeled by values of the pair {AB) are shown in Fig. El The 
operator Z_ consists of the terms 

const x/ §1[ p:V^-HO)I^'Hi.)I^^H^.) .7T_tl . :Cti. - A,B,C^±. (5.27) 
Jc, 27ri Jc 27ri smh smh " 

The integration contours for these terms labeled by values of the triples (ABC) are presented in Fig. |21 Note 
that the contour for the integration variable ^2 is always below the contour for ^1 , which makes it possible to 
avoid the pole at the point £,2 = £,1 + itt. It is interesting to note that for the products S'+S'+y and VS-S- 
the contour for £2 is above the contour for ^1 in the case B — C = — only, where there is no such pole. 
This observation makes these operator products, which enter Eq. (|5.2U|) . well defined. 
Consider now the free particle point p = 2. At this point the action takes the form 

J \ Sir TT Stt / 

where the term — M^x^/Stt in the Lagrangian appears due to the quantum corrections, which make the action 
to be well defined. Hence, it is possible to compare the result with the known results for the sine-Gordon 
model at the free fermion point. 

At the point p = 2 we can derive from the S matrix (|5.5() that 

Z±i0,)Z±{02) ^ -Z±i92)Z±i0,), Z+i9i)Z^{92) = -Z^{92)Z+i0,), 

Zoi9,)Zo{02) = Zoi02)Zo{0i), Z±{0,)Zoi92) = Zo{02)Z±{0,). 



17 




(- + +) (- + -) (--+) ( — ) 




Fig. 3. Integration contours in the terms l|5.27^ consisting the operator Z~{9). The pictures are labeled 
by values of the triple (ABC). The bold points (•) correspond to poles in the variable ^i, while the circles 
(o) correspond to poles in ^2- A circle with a dot in the center (s) means a pole for both variables. 



It means that the operators Z±{9) describe free fermions while the operator Zq{9) describes a free boson. 

Despite of simplicity of commutation relations, the limit p — > 2 in the free field construction is rather 
complicated. First of all, in this limit the operators 

= le'^e^i'^/^):, /(±)(0) = :e-'^(«±'-/2): (5.30) 

are all mutually anticommuting. This expression determines the limit for Z^{9). The situation with Zq{9) 
and Z_ (9) is more complicated. 
We start with the operator 

CE(r, 9U/2 r rif: „fc(e-C)-AK+BK 

{P~l)' sr^ , ,.a+Bt(AB) t(AB) _ / «4 ,,(A)^a\T(B)r^^^'^ 



Zo{k,^\9) c _1^(p2L_ y: (-i)^+^J(^^\ J(^^) = / ^^v^^Ho)i^^Ho 



sinh ■ 



A,B=± 

In the limit p ^ 2 the integration contour in J^"' ^ is pinched between the poles at ^ = 9 + in and 
^ = 9 + i7r(p — 1). The contour for the product in is pinched between the poles at ^ — 9 — in and 

^ 9 — i7r(p — 1). This results in the following identities 

^ e-''^'=-2«.^(+)(5()j(-)/g, ^ j^v ^ f^^^-^j, integral, 

p-2 

^ e-'''''+^'':V'^-\9)I^+'>(9 - in): + finite integral. 

p-2 

The sum of these two terms is equal to 

p — 2 \ J 

It follows from Eqs. (1223 that iTrfc + 2k - (p- 2)^/2 and from Eq. (E^H) that ^ (P - 2)^/^- Hence, 

the sum J^^ ^ + J*^ is proportional to (p — 2)^^/^, which is much greater than the finite contributions of 
J(++) and J(— ). For this reason, we introduce the operators 



bft)-i|^y ^ ' h{t), [b{t)Mt')]^2sinhnt5{t + t'). (5.31) 

Finally, we obtain that the operator Zq{9) depends on the parameter b defined in Eq. H5.23f) only and is 
given by 

/•OO 

Zo{b\9)= dtb{ty''* + V2^b. (5.32) 
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Comparing this with (|5.8|) and (|5.1U|) we see that this correctly reproduces the form factors of the operator e^^ . 

Now we proceed with the operator Z- (9) . The integration contours are again pinched between the poles 
at = 9 zt m and = 9 zt i7r(p — 1) (i = 1,2). This pinching leads to taking out one integration. The 
remaining integration contour goes above or below the point 9 = depending on the kind of the remaining 
screening operator, 5*+ or S-. In the limit p — > 2 the sum of all terms reduces to a residue at ^ = 9. The 
answer for the fermion subsystem turns out to be 

Z+{a\9) = pce'^«(:e'^(''+'^/2).giW2 _i.ei0(e-W2).e-i-/2)^ 

Z^{a\9) = -pc-ie^-e-'^«(:e-''^(«-'-/2):e'-'^/2 +i:e-''^(«+'^/2).g-iW2) ^^'^^^ 

with c = 2-1/277-^6'^^ for p = 2. 

This expression can be compared to the known results for the free fermion point of the sine-Gordon 
model. Though it differs from the known representation [9] for the vertex operators of the sine-Gordon 
model at the free fermion point, it leads to the correct form factors. Indeed, the straightforward calculation 
of the trace gives 

((Z+(z9„) . . . Z+{di)Z^{9n) ■ ■ ■ Z^{9i))) 

a,sinh^^^^M(£^ • ^ - ^ 

It can be proved that this expression is equal to the known one [26,27] 

((Z+(z?„)...Z+(^i)Z_(e„)...Z_(ei))) = (-)"("+i)/2e''S.(^>-«.)(isin7rar:^^ ' ' 



J cosh %^ 



(5.35) 



6. Cosine-cosine model: pi + p2 = 2 case 



The limit pi + p2 ^ 2 {p^ — s- 0) is not singular for the operators corresponding to fundamental particles. 
Nevertheless, the first bound state splits again into a triplet and a singlet. 
Let 

a{t)=a2{t), b{t)=ai{t), p = P2- 
Without loss of generality we shall assume that 

1<P<2. (6.1) 

In the limit ^ they satisfy the commutation relations H5.13|l . The corner Hamiltonian is given by 
(|5.22a|l . As in the last section, we denote by (j) and tp the fields associated with the oscillators a{t) and b{t) 
according to Eqs. H3.21|l respectively. Then 

V3{9) = :exp(i(/)(0; -p) + i^,{9; 2 - p)):, (6.2a) 
(9) = :exp (i^(0; 2 - p) ± i(0(±) {9- p - 2) + ^^^^ {9-p))):, (6.2b) 
4^\9) = :cxp {icj,{9;p) ± i(0(±)(0; 2 ~ p) ~ ^(±)(0; -p))):. (6.2c) 

The fundamental particles are given by Eqs. Ij3.35bll^.35e|l with the operators V3, /f^' given by Eq. (|6.2|l . 
According to Eq. (|2.8|) the mass of the fundamental particles is related with the parameter in the action 
as 

"'=il^= ■ x(2-p) - (6-3) 
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Beside the fundamental particles there is a set of bound states with the masses 



7T(2 — p)r, 

sm — — — 



Mn = 2^i n=l,2,..., n(2-p)<l, (6.4) 

sin ^ ^ " 



which form quadruplets as in general case. 

Consider the quadruplet of the mass Mi = 2/_t. It is instructive to consider the bound state vertex 
operators Bl'^ {i/ = +, S, A, — ) defined in Eq. (|3.8() for small but nonzero p^. In this case we have 



ru-^«»in- "t'-^J^, ^'1' HP3l"y-;^. (6.5) 

The operators i3^'^(6'), Bg'^{d), and B^^{6) correspond to the triplet. The coefficients at these operators 
in Eq. (|3.8|l are finite. The operator B]^^{9) corresponds to the singlet. It describes the field ipsix) that 
decouples and becomes a free massive field in the limit p3 — > 0. The operator B^^{0) explicitly commutes 
with the operators Z^^i {6) . It means that the corresponding particle does not interact with the fundamental 
ones. Besides, it is not a bound state of the fundamental particles, because K^^^ ^ as ^ 0. 

It is remarkable that the operators for the first bound state particle can be written in terms of the vertex 
operators used for the sausage model (see Appendix |E|) , if we assume that 



a{t)=a2{t), b{t)^ai{t), ^/2 cos ^ ^ -i J2| cos || (6.6) 



in the definition of the triplet operators (|5.22fl and 



, sinhTTt /2sinh^^%^ ,^ 
«W = -^y^^a3(0 (6.7) 

in the definition of the singlet operator H5.8|l . Namely, 

B'^:\9) = ic^'Z+{k2,n2\0), bI'\9) = Zo(fc2,K2|0), B^_:\9) = -ic2Z_(fc2, k2|0). 



B'/{e) = z'o{a3\e), 



(6.8) 

where the constant C2 is defined by Eq. (|3.36|l . The scattering matrix of the triplet is the analytical contin- 
uation of the sausage S matrix (|5.5|) to the region p < 2. 

The fundamental particles, the first triplet bound state and the higher (n > 2) quadruplets form the 
particle contents of the cosine-cosine model with the action^ 

>Scc[yi,H d'x (^ ^^^'^^^' + ^^^^^^' +^cos«iy>iCosa2y>2) , aj + al^l. (6.9) 

Note that the form factors for the first bound state of this model can be analytically continued to the 
sausage model region p > 2 by taking a2 — a, ai ^ — i/3. The situation is very similar to that of the 
sine-Gordon and sinh-Gordon model, where the form factors in the sinh-Gordon model are the analytic 
continuation of the form factors of the sine-Gordon model that contain the first bound state particles only. 

Note, that the S matrix of the fundamental quadruplet ^£^£2(0) and the triplet Bl'^{9) {v = +, S, — ) is 
rather simple. Namely, let 

Z,A0i)By{02)= E S^'H0i-02)lCBli\e2)Z,,,{0i). (6.10) 



''This is a corrected (integrable) version of the model proposed in Ref. [18]. 
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Then, the nonzero matrix elements of the S'^^-'(^) matrix are given by 

cosh — , , cosh — 



cosh«±i^' ' ' COSh^^3W2' 

p p 

cosh^ ^''^/^ 
p p 

I J+- V ;+s V w-s V ' p' cosh cosh ^±ii^ 

Consider now the case pi = P2 = 1- At this point the action can be rewritten as a sum of two sine-Gordon 
models at the free fermion point: 



V2 



This free fermion point provides one more check of the construction. In the limit pi —^ 1 the integration 
contours in Zg,^/ (9) are pinched and it is possible to take all integrals explicitly. Let 

a±(t) =al(^)e±"l*l/4±a2(^)e^"l*l/^ [a±{t),a±{t')] ^ tS{t + t'), [a+(t), a_(<)] = 0. (6.12) 

Define 

ct>±{0)= / ya±(i)e'^*. (6.13) 



Then 



where 



Z++{ki,k2\0) = uje''+^ re'-^+^^^i, (6.14a) 

Z+_(fci, fc2|0) = iope'^^/^e''-' (e'^'=-/2 le'^-^^+^^i - ie-'^'^-/^ le'^-^''-^):) , (6.14c) 

Z__(fci, fc2|6') = -tje^^e-'=+^ (e-''^'=+ :e-''>+(«+'^): + :e-'"^+(^-''^):) , (6.14d) 

fci ± fc2 ai ± a2 , , 

fc± - ^ - (6-15) 



and the operators give the form factors of the exponential fields 

Oaiaiix) = e'°i'^i(^)+''*2V2(2;) _ gifc+xi X2 (a;) ^ (6.16) 

The vertex operators anticommute: 

Ze,e[{Ol)Z,^e'S02) ^ -Z,^,,i02)Z,^,,i9i), for pi = p2 = 1- (6.17) 

The pair Z^+{6), Z [9) describes the Dirac fermion related to the field xi, while the pair Z |-(6'), Z^ [9) 

describes that related to X2- The representation (|6.14a|l . (|6.14d|) coincides with that described in [9] for the 
field e""^ (see Eq. (|5.28l) for normalizations) with a = k+, while the representation (|6.14b|) . (|6.14c|) gives the 
representation 15.33|l with a — k-. Both of them provide the right hand side of Eq. H5.35|l with a = fc± as 
the final result for form factors. 
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7. Conclusion 

The form factors of the exponential operators (|3.11() in the model with the action H2.1() in the unitary regime 
can be expressed as traces of the vertex operators, which are realized in terms of three sets of free boson 
operators. This bosonization procedure provides an algorithm for construction of an integral representation 
for form factors. For the fundamental particles a general explicit integral formula (|3.15l) - (|3.19|l for any 
iV-particle form factor has been found. 

There are two important limiting cases. The first one, p2 — + 0, gives the sausage model in the dual 
representation. The free field representation for the vertex operators is constructed by the limit in the 
free field representation for the first bound state of two fundamental particles. The second limiting case. 
Pi + P2 2 gives the known cosine-cosine model, that contains the fundamental particles as well as their 
bound states. The first bound state forms a triplet and the respective vertex operators can be considered as 
analytic continuation of the vertex operators of the sausage model after the substitution pi P2, P2 —^Ps, 
P3 Pi- 

There are three interesting problems, which we did not solve in this paper. The first one is to compute 
analytically the integrals for two-particle form factors of some special exponential operators. By the analogy 
to the sinc-Gordon theory, it is expected to be possible to compute them for Ui — {n € Z). In this paper 
the integral was only computed for the one-particle form factor in the sausage model (see Appendix IB|| . 

Another task is to obtain a free field realization of the vertex operators for restricted models like the 
parafermion sine-Gordon or the coset models, which are related the family of models under consideration 
by the quantum group reduction. It is expected that it can be achieved by an appropriate change of the 
screening operators [28] together with, probably, some Felder type resolution. 

The third and the most important problem is related with the prescription for calculation of the form 
factors of the descendent operators H3.10|l and of the nonlocal with respect to the fields (fi{x) operators. 
These problems will be addressed in forthcoming papers. 
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Appendix A. Regularized t-integration 

The Jimbo-Konno-Miwa regularization rule, described in Sec. El leads to the following simple formulas for 
regularized integrals: 




(A.l) 



(A.2) 



(A.3) 



(A.4) 



Here Ce is the Euler constant. 
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Appendix B. List of trace functions 

Here we list formulas for the constants Cu^ and functions GjjjUki^) defined in Eqs. H3.39a() . Ij3.39b|l . which 

determine the traces according to Eq. (|3.38() . for the operators Vi{9), l[^\9) defined in Eqs. (|3.23() and 
V'^^\e), l'^^\e) defined in Eqs. ^T^. These calculations are performed straightforwardly by application 
of the formula (we recall that subscripts «, j are defined modulo 3): 

sinh^ — 1 
«°''"°^"'»-' s.„l.,i...hM^l-e-^.. ^"^" + '''- 



For the constants we have 



dt e-"* sinh^ ^ sinh "^+1+^'+^)* 



Cy.=C^ = --^[-l fssinh^^sinh^sinh^j' ^^''^ 

-1/2 

^4*'"^*" r(i + f/p,)' ^ ' 

Cy{±] = Cj(±) — C — Ci\p^^p. (B.3) 
For the functions Gjk the result is 

Gv^vM^G,,{e), (B.4) 

Gy,^,±m = G^,^>^ie) = wl^He), (B.5) 

GjiA>jiB,{e)=G\f\9), (B.6) 

Gv(A,v(B) W = G,(^„(B,(0) - &^^He) = G^f^H0)|p,^p, (B.7) 
Gy^AUiB,{0) = Gi,A,y,Bm - = {&^^\e))~\ (B.8) 

The functions Gij [9) are rather complicated: 

( r°° dt sinh^ ^ sinh ^ A 

= =p (- r '1 "-"j;-'^:-" ) (. ^ . M ...). (B.9b) 

\ Jo t smh 7rt smh / 

but they never appear in the integrands, i. e. their arguments never contain integration variables. The 
functions W^^/' (0) can be expressed in terms of a single function W{p; 9): 

W^f\e)^l, (B.lOa) 

w^%{d) = w{p,-e±-ifp,), (B.iob) 
w'>^l{e)^w{p.c,e), (B.ioc) 

, ^ dt sinh ^ 003(6* + m)t 

The function W^(p; 6*) possesses poles at the points = if{p - 1 + 2Mp) + 27riiV (M, TV = 0, 1, 2, . . .) and 
6* = - fip - 1 + 2Afp) - 27ri(iV + 1), and satisfies the relations 

W{p;0) = VF(p;-6l- 27ri), 
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The functions g[^^^ (0) are purely trigonometric: 

G^^^\0) ^ ±2isinh^-i^, (B.12a) 

Pi 

Gf^\e) = 2tanh-sinh^-i^, (B.12b) 

2 Pi 

G%tli(^)^G%tl{6)^l, (B.12C) 
G^iO) ^ G^lie) ^ -icoth ± . (B.12d) 

It means that the form factors of the sausage model can be expressed in terms of multiple integrals of 
trigonometric functions with two different periods. 

It is important to make the following remark about the contours. To fix the contours it is necessary to 
consider the product in the trace and look at the poles of this expression. There are poles that appear from 
the operator products (|3.24|) . H5.15|l and the extra factors in the screening operators. The prescription for 
the integration contours with respect to those poles was already described in the main part of the paper. 
The additional poles that arise due to the trace should be treated as follows. Namely, the contours should 
go below the excessive poles at the points Oj + \5j with 5j > and above the excessive poles at the points 
9j + i6j with Sj < 0. For example, the trace 



sinh sinh 2 sinh sinh 



p p p p 

possesses poles at ^ = 6* ± Itt, though the product 



np 



possesses a pole at ^ — 6 — in only. It means that the integration contour for ^ must go below both poles at 
9 ± m. In contrast, in the expression 



sinh sinh 2 sinh sinh 

p p p p 

the integration contour must go above both poles at ^ = ± Itt, because the only pole at ^ = + Itt coming 
from the extra factors survives after removing the trace sign. 

As an example of application of the integral representation, we calculate the one-particle form factor 
{{Zq{9))) in the sausage model. This form factor 

_ {Q\e''^^+'^\9,Q) 
- (ole-^+fxlo) 

is just a constant. Nevertheless, it gives the first nontrivial contribution to the infrared asymptotics of 
correlation functions. 

Using the fact that W^++'>{9) = W^-~^{9), we obtain 



{{Zo{9))) = -0^2 cos ^ J ^ ((t/We'« - iy(-)e^) (/^e'' - i/^-^e"'')) - 



sinhS^^ 

p 



p \J._ 27ri sinh 7c-+ 27ri sinh ^ 



The contours and C |- are chosen according to the above mentioned rule. As we have seen just now, 

the contour C+_ goes above the poles at ^ = 9 ±m, while the contour C |- goes below the poles at these 

points. 
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With the substitution ^ ^ ^ + 6* we obtain 



2tt sinh ^ sinh £^ 



The integration can be done by residues, if we suppose e. g. Imfc < 0. With these integrands, which are 
periodic functions times an exponential function, it is easy to sum up the residues with the result 

„l/2 sin ZL si^TT^ 

((Zn(fl))) ^ P - s^simrPb, p^2a^ = 2ff + 2. (B.13) 

sin — TTaa 
p 

In the limit p ^ 2 we have ((Zq -\/27r6 in consistency with the free boson result. Note, that the 
result ljB.13|) is invariant with respect to the reflection transformation 6 — > 1//3 — b (see Ref. [16]). 

Appendix C. Check of the Zamolodchikov— Faddeev algebra relations 

The relation (|3.1c|l for the commutator [H, Z^^i (9)] of the corner Hamiltonian and a vertex operator can 
be checked straightforwardly. Here we verify the commutation relations H3.1a|l of two vertex operators and 
prove the expression H3.36|l for the normalization constants by checking the relation IjS.lbp . 
We start with the simplest commutation relation 

where 9i2 ~ 9i — 92 and we used the fact that the S matrix Sp-^p^{9) is a tensor product of two sine-Gordon 
S matrices (EUDl)- Substituting Eq. I|3.35b|l we obtain 

^3(01)^^3(^2) = -'Sp,{9,2)llSp,{Sl2)llV3(92)V3{9i). 

Reduce the l.h.s. and the r.h.s. to normal products with the help of Eq. H3.24a|l : 

533(-ei2) :^3(0l)V^5(^^2): = -^Pi(^?i2)it^P.(0i2)ii533(f?i2) :^3(0i)F3(02): 

i.e. 

-Spi{6)XtSp2{S)t+ ^ '^^^ in\ ■ (C-1) 
333 (t/j 

The last identity can be obtained using Eq. H2.10|l . 
Consider now another commutation relation 

Z++{9i)Z^+{92) = -Sp^{9i2)X_Sp2{9i2)\\Z^+{92)Z++{9i) - Sp^{9x2)j^tSp2{9i2)\\Z++{92)Z^+{9i). 

Using Eqs. JOHll . (TOTll . JiHlall . i^Hbll . (TOl we obtain that it is equivalent to the following equality 

.i,,a ,fiA)r.,.93,{ei2)wii\^-9,)wii\^-92) 



sinhi^^^^ 

pi 



^ SpA0x2nzspA9,2nt E / ^ Mo.)%{e2)i[^\o: ^33(-^^^)-^i;(^^ - - 



SpA0i2ntSpM2)XXT. J §r-Moi)%(02)ii^\o 



with 



4-6)i-ij/2 



sinh 

pi 



V,i9) = y,(0)e^^^4^^ /W(C) = /^^(O = -ilt^iOe--' 
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We note that operator part equal to '■V3{0i)V3{92)i\^\^): is the same m all terms, therefore it is sufficient 
to prove the identity for the remaining functions. Here and later we need the identities 

J^%i-0) sinh^ 4tli-9) _ cosh^ _ sinh^ ^^^^^ 



4%(0) sinh^' cosh^' g^-\0) sinh^- 

Dividing the functions in the l.h.s. and r.h.s. by the numerator of the function in the l.h.s. and using Eqs. H2.9|l . 
(|2.10|l and l|C.l(l we obtain 

1 sinh(j/i — 2/2) sinh(x — yi + z/2) 1 



sinh(a; — j/2 ^ -^Z^) sinh(z — yi + j/2) sinh(a; — j/i — z/2) sinh(a:: — 1/2 ~ -2^/2) 

sinh z 1 



sinh(a: — yi + j/2) sinh(a; — yi — z/2) 

with X = ^/p, yi — 9i/p, z = in /p. This identity can be checked straightforwardly. 

The other commutation relations can be checked in the same way. The commutation relation 

Z-+i6i)Z++{62) = -Spj^{9i2)l+Sp2{6i2)X+Z-+i(^2)Z++{9i) - Sp^{9i2)_XSp2{6i2)X+Z++{(^2)Z-+{9i) 

is valid due to the identity 

1 sinh(a; — 2/2 + 2/2) sinhz 1 sinh(a; — yi + z/2) 

sinh(a; — yi — z/2) sinh(a; — y2 — z/2) sinh(z — 2/1+2/2) sinh(a; ~ y2 — z/2) sinh(a; — yi — z/2) 

_^ sinh(yi - 2/2) 1 

sinh(a; -2/1+2/2) sinh(a; - yi- z/2)' 

The commutation relation 

Z_+(0i)Z_+(02) = ~SpA0i2)ZZSpM2)XtZ-+{02)Z-+i9i) (C.3) 

contains a two-fold integral. There is an important fact concerning the integration contours. Each product of 
vertex operators in this equation splits into four terms containing different operator products V{9i)V{92) x 
I^^K^i)I^^KC2), A, B = ±. For each pair {A, B) we can chose the contour cj"^'^-* for ^1 and the contour 
^(-4,5) jT^j. that it satisfy the pole avoiding rules for all three operator products in Eq. ljC.3|) . The 

important fact is that c[^'^^ can be deformed into €2^'"^^ without meeting poles. (For A = i? it simply means 
that C^^'^^ and C^^'^^ are up to a deformation the same contour.) This means that we can symmetrize the 
integral in the integration variables. After this symmetrization we arrive to the identity 

V sinh(.X2 — xi — z) , 

f{xi,X2)+f{x2,xi)—^ ^—4=0 (C-4) 

sinh(.T2 — xi + z) 



with 

/(a;i,a;2) 



1 sinh(a;i — 2/2 + z/2) 



sinh(a;i — i/i — z/2) sinh(a::2 — 2/2 — 2/2) sinh(a::i — y2 — z/2) 
1 1 sinh(a;i — 2/1 + z/2) 



sinh(a;i — 2/2 — z/2) sinh(a;2 — yi — z/2) sinh(j:i — 2/1 — z/2) ' 

which is checked straightforwardly. 

The relations proved above contain the vertex operators Z^^' (9) with e' = + only. That is why the 
proofs only involved the screening operator Si{9). The check of relations that contain the vertex operators 
with £ = + go essentially in the same line. The proofs of these relations have to do with the screening 
operator S2{9) only, and the identities to be checked are obtained from the above ones by the substitution 
yi — * yi + i7r/2. The remaining relations for the vertex operators with generic e, e' follow from the proved 
ones and commutativity of the screening operators Si(9) and 82(9) . 
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At last, we check the normahzation condition (|3.1b|l . It is enough to verify it for the product 

A.B—:L ^1 ^ C2 



As a rule, the pole at f?' = 9 + 'nr is known to be related to pinching contours between poles [7]. An elaborate 
analysis of the poles shows that the only pinching situation takes place at the term with A = B = +. The 
contour Ci goes below the pole at 6*' — ^ coming from ^3^' and above the pole at ^2 + coming from 
9it^'' ■ The contour C2 for ^2 goes, in turn, above the pole at 9+^ — from Trwgf ^ (^2 — ^)/ cosh ^^"^^'^ ■ 
We have to push the contour Ci, e.g., upward. The residue at ^1 ^ 9' — ^ contribute the singular term, while 
the remaining integral along the contour that goes far from poles does not contribute it, because ^2 will not 
be pinched between the points ^1 — and 9 + ^ — i^^. In the residue we push C2, e.g., downward. The 
residue at ^2 = ^ + 3- — contributes the singularity. 
For the operator part it is easy to get 



:V{9 + m)V{9)i[+\9 + if )/^+' + f - = 1- 



Finally, we obtain 



( pi ) ^ ( P2 ) 



9 -9-m r (1 + ^ ) r (1 i > 



pi 7 V P2 



where Vp is defined in Eq. (|3.26d|) . Comparing with the normalization condition Ij3.1b|) we obtain the 
value of the product C1C2. This calculation does not fix ci and C2 separately. But if we demand all the 
operators Z^^,{9) to be normalized by Eq. HB.lbjl . we get all the products QCi+i (which are the result of 
cyclic permutations of subscripts in the expression for C1C2) and obtain Eq. 



Appendix D. Vertex operators in the sausage model 

Here the operators Z+{9) and Zq{9) will be obtained directly from the general model H2.1|) by taking the 
limit P2 ^ 0- The operator Z^{9) will be found from the commutation relations of the ZF algebra for the 
sausage model. The proof of the commutation relations in this case has some additional complication related 
to the pole at ^ = in the product V^^\9)I^^\£^). It will be shown how to manage this pole. At the end, 
the normalization constant c in Eq. H5.22|l will be derived. 

Consider the general model H2.1II for P2 < 1 and calculate the products 

Z,+ {9' + if (1 -p2))^e'-(^ - f (1 -P2)) (D.l) 

in the vicinity of the point 9' = 9. The residue of this pole will give the first bound state of two fundamental 
particles. 

We begin with the product 

Z++i9' + 6)Z+_{9 - S) 



-PC2 



_± V,[9' + 5)V,{9 - S){I^-^\Oe^^ - i/ "^(Oe"''^) — 

27ri smh ^ 



P2 



-pc2y g.:y3(e' + '5)F3(^?-^)(/W(e)e- -i/^)(e)e--):e^(«+^')-^2« 
X g,s{9 ^9'- 2S)wit\^ -9'- 5)wit\^ - + S)- 



sinhi^ 

P2 
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with (5 = ~P2)- The integration contour is pinched between the pole at the point ^ = 9' oi the function 
'w{p2^ 1/2|^ — 9' — 5) and that at ^ = of the function sinh""'^ Pushing the contour downward and 

calculating the residue of the second pole we obtain 

Z++{9' + 5)Z+^{9 ~5) = ■.V:,{9 + 5)^3 - <5)(4+^(0)e'== - ui-\9)e-^-):e^^' + 0(1), 

^ ^ \/7rG'(p2,-i7r) 

G(pi,-2,5)G(p2,-2<5)- 

It is straightforward to find that 

■■Vz{0 + .5)^3 - 5)1^2'^ {Oy. - V^+\9l 
■.Vz{e + 5)V:i{9 - 5)i-\9): ^ 

as (5 — > i7r/2 (p2 — *■ 0). In this limit the constant C remains finite: 



Here we used the identity 

As a result, we obtain 
It is easy to check that 



Gfe-ivr) _ 2 r(l-f) 



(D.2) 



G{p;-iTr + mp) ^7:{l-p) T {^) 



Z++{9' + 5)Z+_i9 -S)^ ^^pVik,, + 0(1). 



c 

— = r2,i, P2 = 0, (D.3) 
c 



with the constant c given by Eq. (|5.25() with p ~ pi and r2,i given by Eq. (|5.7(l . It proves that -B_|_' {&) 
c^^Z+{k,K\9) with Z+{k,K\9) defined by Eq. (|5.22b|) . 

Consider now the products ljD.l|) for e' + e = 0. For e' = +, e = — we have 



z++{9' + s)Z--{e-s) 

_^2g^^ii±a(e+9')-fei5i-fc252 



ciC2l^3(e' + <5)^3(f?-^) / ^ / ^(/l+^(ei)e''^-i/^^(ei)e-''^)(4+^(6)e^^-i4"^(6)e-''^) 



sinh ^^''-'"P^/' sinh Si 



ClC2 



i: i Si S:*3(»' + fl»3(»-.)/!-fe)/f'te): 

13 y sinh sinh 

Pi P2 

There are two pinching points as 9' 9. The first one is the same as in the case e' = e = +: the contour 
Ci in this case is also pinched between the poles at the points 9' and 9. The second one is in the term 
A = B = — and involves both integration variables: the contour Gi goes below the point 9' — and above 
the point ^ + while the contour G2 goes above the point 9 — mp2- We want to move the contour Gi 
so that it would go above the points 9 + i7rp2/2 and 9' — mp2/2. So we have to push it through the pole 
9 + mp2/1 in the term with A = +, i? = — , and through the pole 9' — mp2/'2, in the term with A = B = +. 
After all these operations we obtain 

Z++{9' + S)Z__{9 -S) = ^^{Jo + Ji + J2 + J3) + 
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with 



Ji = -C'^:V{0 + 6)V{e - S)l[+\0 - 

J2 = -c'-^:V{e + s)v{e - 6)ii'^\e - '^)ii-\e):, 
rO-Tf)r(i- 



J3 = 



-2 \ Pi J \ Pi 



where the integration contour in Jq goes below the point 6 — in + inp2/2 and above the points 6 + i7rp2/2, 
6 — i7rp2/2. Similarly, for e' = —,£ = + we have 

z_+{e' + s)z+_{e -s) = ^^{4 + J( + 4 + 4) + o{i). 

Here the operators Jq and J^' (i = 1, 2, 3) have the form 



4 



■ rll^ Sinh g-^+i^(2-P2)/2 



pi 

sinhMi^ 



J'^J, El (i = 1,2,3), 

sinhi^L ^ ' ' ^' 

pi 



where the integration contour in Jq goes above the points 9' + in — i7rp2/2, 9 + i7rp2/2, 6* — i7rp2/2. 
Now we should take the limit P2 — » 0. In this limit the constant C remains finite: 



^/ -2 V PI / n 

It means that 

Jl + J3 = V^C"piP2"'^'4(^) + 0(P2^'). 

We want to obtain an expression for B^^{6) in the limit P2 — > 0. We should consider the difference 

Z++{9' + S)Z__{e -S)- Z_+{9' + 5)Z+_{9 - 5). 

Due to the relation 

sinh ^-«+'n2-i^.)/2 2sinh^cosh^^^|^ 2 sinh cosh 

sinh^-^-'"(^-^-)/^ sinh^-^-'"<^-^-)/^ sinh«^ ^^^'^ 

pi pi pi 

we have 

ciC(Jo - J'o) = - — ,/2cos— Zo(0) + 0(p2) = -rifi'iZo(0) + 0(p2). (D.4) 
c V f^i 

Similarly, from the relation 

sinh^^fasl 2 sinh ^ cosh ^^^1^:^ 
1 ?^ — = ^ = 0{P2) 

smYi'-^ sinhi^ ^ ^ 

pi pi 
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we obtain 

Hence, the only finite contribution comes from the difference Jq — Jq. The equations ljD.4|) . IjD.Sp prove that 
B^/{e) = -Zo{k,K\0) with ZQ{k,K\d) defined in Eq. f^?^ . 

To obtain an expression for Bg'^{9) we should consider the sum 

z++{e' + 5)z__{e -5) + z_+{e' + 5)z+_{e - 6). 

It is easy to check that 

sinh 2coshS^sinh^I^^|^ 2 cosh sinh i^L 

^.^^ .^-g-i.(2-p.)/2 3.^t^ C-^~i^(2-p.)/2 Sinhi^ ^ ^^'^ 

Pi Pi Pi 

and hence, 

Jo + = 0(1). 

Similarly, from the relation 

sinhi^iizM 2coshi^sinhi^^|^ 

1 + = ^ = 2 + 0{p2) 

sinh^ sinh^ ^"^'^ 

we have 

Ji + Jz + J[ + J3 = 2V^C"piP2"^^^^o(^) + 0{p\'^), J2 + J2 = 0(1). 
It is straightforward to find that 

1™ ^.rz„r<'„. — 



lim 2\/ttcC' pi ~ lim py^i^o'^ 

P2^0 P2^0^ ^ 



Together with Eq. I|D.3|I it proves the expression (|5.8(l for i?g'^((?). 

Derivation of the operator Z- [9) by taking the limit of Eq. I|D.1() for e' = e = — would be very compli- 
cated. In fact, we do not need it at all. It is natural to suppose that the operator Z^{9) has the form of a 
linear combination: 

c++sl{e)v{e) + c+^s+{e)v{e)s^{9) + c^+s^{9 - 2TT\)V{e)s+{e + 2tt\) + c^^v{e)sl{9). (d.5) 

All four terms in this combination are different: though the integrands are identical, the integration contours 
are not the same. To see it, we consider, for example, the difference 

A(6i) = S+{9)V{9)S-{9) - Sl{9)V{9). 

We can expand it as 

A(0) = A<^^^^H9), 

A,B,C 

where 



Jc^ 27ri Jc_ 27ri smh smh ■ 



Jc+ 27ri 7c+ 27ri smh ^ 



V(+) (9) = iy W (6»)e-''+'=^ !/(-) {9) = F^-' (e*)* 



p 



sinh ■ 



The operator A'^^^'^'^9) is a difference of two terms, which are distinguished by the integration contour for 
the variable ^2- For this reason the calculation amounts to finding the residue at the point ^2 = 0: 

^' 6=0 ic^ 27ri ^ sinh Si^^ sinh Sa^^ 
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The only nonzero contribution comes from the terms with A = C. As in this case 

Resl/(+)(0)/(+)(O = -Resl/(-)(0)/(-)(O npr-'^-}^, (D.6) 

we obtain that integrands for the operators A^+++)(0) and — A*^ ^ ^(0) coincide. The same fact takes place 
for the operators \9) and — A'^"' If we consider the contour in the variable for the last pair, 

it is the same for the operators A^ )(6') and A^+^+\e). It means that A( ^9) + M+-+^{e) = 0. We 

have a different situation for the operators A(+++'(6') and A^ ^ H^). In this case, both integrands possess 
poles at ^1 = 9 — Itt, but the origin of these poles is not the same. In the term A^ ' ^ (9) it comes from 
the product V^^\9)I^^\^i) and the contour for the variable must go below it (see Fig. In the term 
A(+++)(6') the pole arises from the product and the contour for the variable must go above 

it. It means that 

A{9) = A(+++)(0) + A(-+-)(e') = Res Res T>("H6')/(+)(?i)/*"H6)- 



6=e-i7rC2=e ^ ' ^"'sinhii^^^ sinh^^^^^ 

p p 

2^ 



sinh ■ 



p 



It proves the first line of Eq. H5.2U|I . The second line is proved similarly. 

To fix the coefficients Cab we have to analyze the commutation relations for Zi{9). It will be shown 
below that the only solution consistent with the commutation relations is the solution (|5.22d|l . i. e. 

C+_ = const , C— ^ C++ ^ C-+ = 0. (D.7) 

Now we can prove the commutation relations of the ZF algebra for the sausage model. This problem 
is more difficult than in the case of the operators Z^^i{9). The proof of the relations for Z+{9i)Z+{92), 
Z+{9i)Zq{92), and Zq{92)Z+{9i) is, in principle, the same and reduces to the identities 

S{0)XX = -93z{e)/g33{-9). 
1 sinh(yi — 1/2) 1 sinh(a; — yi + z) 



sinh(a; — 1/2 — z) sinh(?;i ~ 1/2 — '2z) sinh(a; — y2 — z) sinh(a; — yi — z) 

sinh 2z 1 



sinh(j/i - t/2 - 2z) sinh(a; - yi ~ z) 

According to Eq. (|5.19|l a single contour can cross the pole a± — 9. Hence, the contour in the 1. h. s. of the 
commutation relation can be deformed into that in the r. h. s. The situation with 

Za{9^)ZM = S{9i - 92)+oZ-{02)Z+{9i) + S{9i - 92)-„ + Z+{92)Z^{9,) + S{9i - 92fo°oZo{02)Zo{9i) (D.8) 

is more complicated. First, as it was mentioned in Sec. |S1 the contours in the definition of the operator 
Z-{9i) cannot cross the pole at 9i. Second, the contour in the 1. h. s. in the operator Zq{9\) goes above the 
point 92 and the contour in Zq(92) goes below 9\. On the other hand, in the r. h. s. a contour must go above 
the point 9\ if it belongs to the operator ^7(^2) and below 6*2 if it belongs to Zi(9\). 

We shall prove this identity in two steps. On the first step we shall consider all integrations so as if they 
would be done along the same contour and show that the symmetrized integrands of the 1. h. s. and the 
r. h. s. coincide. On the second step we shall move all contours to a single one and show that the residues 
due to the encountered poles cancel each other. 

For any product X of the operators Zi{9\) and ^7(6*2) denote by \X\ this product, but with the contours 
in the screening operators going for example above 9\ and below ^2- We denote as 

A{X\=X-\X\. 

On the first step we prove that 

\Z^{0i)Z^{Q2)\ = S{9, - 92)+^[Z^{92)Z+{9,)] + S{9i - ^2)0^ [^+(^2)^- (^i)] 

+ S{9i-92foo[M02)Zo(9i)]. (D.9) 
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On the second step we prove the equation 

+ SiOi - d2fo°oMZo{92)Zo{ei)]. (D.IO) 



Combming Eqs. (|D.9|) and (jD.lOp we obtain Eq. IjP.Sp . 

The first step is performed just as in Appendix lO and reduces the proof of Eq. (|D.9|) to a check of the 
identity of the form IIC.4|l with 

-2coshz sinh(a;i - ?;2 + ^) 
f{xi,X2) = ^-r- -^—7 T^TT T 

smh(a;i - 2/2 - z) smh(a::2 - y2 - z) smh(xi - y2 - z) 

sinh(2/i — 2/2) sinh 2z 1 
sinh(j/i -y2- z) sinh(?;i - 2/2 - 22) sinh(a;i - 2/2 - 2) sinh(x2 - y2 - z) 
sinh(a::i —2/2 + 2) sinh(x2 —2/1 + 2) 
sinh(a;i — 2/2 — 2) sinh(a;2 — yi — z) 

sinh(j/i — 2/2) sinh 22 1 



sinh(2/i — 2/2 — 2) sinh(2/i — 2/2 — 22) sinh(a;i — 2/1 — 2) sinh(x2 — 2/1 — 2) 
sinh 2 sinh 22 + sinh(2/i — 2/2) sinh(2/i — 2/2 — 2) —2 cosh 2 

sinh(a;i — 2/1 — 2) sinh(.T2 — 2/1 — 2) sinh(a:i — 2/2 — 2) sinh(a;2 — yi — z) 

sinh(.xi -2/1+2) 

sinh(a::i — 2/1 — 2) 

This computation would lead to the proof for operators with the same contours. In our case, however, 
the contours are different and on the second step we should take this fact into account. 

Now we calculate A[Zo(0i)Zo(02)]- We have to proceed in the same way as in the calculation of the 
operator A(6') above. We skip the complete analysis and write the contributions that come from the defor- 
mation of the contours. We have to push the contour Ci attached to Zq{9i) down across the point 6*2 and 
the contour C2 attached to ^0(^2) up across 9i. The result will have the form 

A[Zoi9i)Zoi92)] - -2c2cos- 

P 



+ Res Res V^-\9i)i^+\(i)V{92)I^~\(2) 7^^^ . ^ 

Using Eq. I|D.6|1 and the commutation relations (j5.18|l we obtain 

TT -/ N - ^ ^, ^ 1 sinh — 

A[Zo{9i)Zo{92)] = ^2Tr^c'Dpcoa - {I^-\92 + i7r)T/(0i) - V{92)I^+H9i - i^))- 



^ " ■ ' sinh ^i^^ sinh ^1^^^' 

Similarly, for the terms of the r.h.s. of Eq. (|D.10|) we find 

A[Zoi92)Zo{9i)] = 0, 

A[Z_(^?2)^+(0i)] = c^ Res Res y(+)(0,)/W(a)/(-)(6)m)--^^^--l^^ 
^1=92 C2=92+i7r smh ^ smh ^ 

p p 

= -n^c^Dpi^-\92 + iTT)V{9i) ^ 



sinh i2: ' 

p 



A[Z+((?2)^-(0i)] = c2 Res. ResV(~H92)i^+\^,)i(-H^2)V{9i)—-T^_ 

f 1 — — ITT f — cmh 



TT 



i-C=e2 sinh '^^-^^-'" sinh '^^"^^"'" 

p p 

7T^C^Dpi^+H9 - i^)l>(02)-^ 

smh — 

p 



32 



for the operator Z^{9) defined by Eq. Il5.22d|l . i. e. for the case Eq. (|D.7|) . Taking into account Eq. 1)5. 5|) 
for the exphcit form of the S matrix of the sausage model, we obtain Eq. (|D.10I) . which completes the proof 
of Eq. HU.8|I . The proof of this identity is a basic one, because it fixes the form of the operator (|5.22d|) 
from that of the operators (|5.22b|) and (|5.22c|) . Indeed, if we would allow all terms in Eq. I|D.5|) the value 
A[Z± {9 2)2^(61)] would contain some extra terms, which are absent in A[Zo(0i)^o(^2)]- 

The commutation relations for Z±{6i)Z^{92) can be proved in the same way. For Zq{9i)Z^{92) and 
Z-{9i)Zo{92) the proof is a little more complicated: it remains one integration in products like A[Zo{9i) x 
Z_(02)]- First, we need to make sure that the integrations in both sides of the equation are taken along 
the same contour. Then we have to check that the integrands in both sides coincide identically. One more 
complication appears in the commutation relation for Z-(9i)Z^{92)- There is a double integration in this 
case and we have to symmetrize the integrands in the integration variables before comparing them. We have 
checked all these relations, but we omit the calculations as they are very cumbersome. 

The normalization constant c in Eq. H5.22|l can be calculated as follows. Consider the operator product 
Z+{9')Z-{9). It can be shown that the only contribution to the residue at the pole 9' = 9 + iir comes from 
the product 

2 f d^2 



2tt 2tt ' ' ^^"^sinhii^^sinhSi^^ 

p p 

I ^^■■vi-He')v(+H9)i(^Hi,)i(-H^2y. 



^ 2tt 2tt 



-9^-+H9-9rg^+-H^2-^i) 



g(~+){^, - 9')g(++){^i - 9)g(—H^i - 9')g(+-){^i ~ 9) sinh ^i^^ sinh ^ 



There are two efi^ects here. First, there is pinching of the contour for the variable ^1 between the points 
^1 = 9 and £,1 = 9' — \tt and of the contour for ^2 between the points ^2 = 9 + m and ^2 = S' . This gives 
a double pole in 9' — 9. Second, the function g'^ ^\9 — 9') has a zero at 9 — 9' = —in. It decreases the 
multiplicity of the pole by one. Finally, we have 

Z+{9')Z_i9) = -c^g(-+)'i-in)i9 -9')- ^ 



' gi-+){9 -9')g(—){9 



X Res Res ^ , ^^^^^ — , , . + 0(1) 

gi^ega^e+iTT g(++)(gi ~ 9)gi+-){i2 - g) sinh sinh 



9' -9-17: P ^ ^' r(l/p) 
Demanding that the coefficient at the pole —i/{9' — 9 — 'nr) was equal to 1, we obtain Eq. H5.25|l . 

Appendix E. First bound state in the cosine-cosine model 

Here we calculate explicitly the vertex operator B^^{9) corresponding to the highest component {Q2 = 2) 
of the quadruplet of mass Mi . It appears in the product 

Z++i9' + if (1 - pi))Z^+{9 -fil- pi)) (E.l) 

in the limit 9' 9. Using Eqs. (|3.24a(l . I|3.24b|l . (|3.26b|l . 13.26d|l . we obtain for this quantity: 

r;^^g33{9-9' -26) I ^ (^.V{9' + S)V{9 - S)li+\0:e^'+^'''+'''^^'+''^/^-'''^ 

^[^-2) ^ f iii-9) ,1 l\^f i(e-g) , V 



( i(C-9') _ 1 I \ npi 2 pi 



npi 2 
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+ i:V{e' + 5)v{e - c-'=i+('=i+'=^)(^+«')/2-fci« 

ri^^ + 2j 3 i\ / i(e-e) 1 



r 



r - 



7rpi ^ 2 ' pi 



with S = ^(1 — pi) and = e(CE+iogTp)/p^ xhe integration contour is pinched between the poles at 
£, = 9' — i7rpi/2 and £, = — 'nrpi/2 in the first term and between the poles at ^ = 9' + i7rpi/2 and 
^ = 9 + mpi/2 in the second one. Pulling the contour through a pole in each term we obtain for the 
quantity (|E.ip : 



-yplc,r^,'933{~2S) \ (■.V{9 + 5)V{9 S)l[+\9 izmi); e-'^3+fe.« 
; r ( ^ 



+ r.V{9 + 5)Vi9 ~ S)l['\9 + iiEiy, e'''+''^^^ + 0(1). 



The first normal ordered product here tends to V^~^{9), while the second one tends to V^~^''{9) defined in 
Eq. (|5.14a|) . Recall that though the expressions (|5.13|l - (|5.17f) were introduces for the sausage model with 
p = Pi > 2, they can be analytically continued to the cosine-cosine model case p — P2 < 2. Finally, in the 
limit P3 — > we obtain 

Z++i9' + 5)Z+_{9 -S)^ .^ipV{k2,n2\9) + 0(1) (E.2) 

with V{k,K\9) defined in Eq. H5.16|l and 

^ ^1/2 G{pi,-m) 

G(pi,-i7r(l -pi))G(p2,-i7r(l -_pi))' 



Using the identity ljD.2|l and 



we obtain 



G(p;i7r-mp) ^ p T (f ) 

C2C" 



~Psm-- \ —— (E.3) 



ri,i, Pi+P2 = 2, (E.4) 



where c is given by Eq. (|5.25|) with p = p2- The expressions (jE.2|) and (|E.4|) prove the representation for B]^^ . 

The proofs of representations for the operators i3g'^ and B]j_^ are in the main features the same as for 
the operators i?^'^ and B^_l^ in the case of sausage model and we will not reproduce them here. 

Appendix F. Regime II 

In this appendix we consider the model in the nommitary regime II, where 

Pi,P2,P3>0, pi+p2+p-i = 2. (F.l) 

In this regime the theory is completely symmetric under the permutations of the pairs {p2,^2)i 
{P3,(p3) and possesses the U{1) x U{1) x U{1) symmetry described by the topological charges 

Q^^jdx^]l ^ ^e^-a.^,. (F.2) 

The values of the charges must satisfy the conditions 

Q^ e Z, Q, + Qj e 2Z. (F.3) 
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We would like to study this regime by means of the free field representation technique. In Sec.^we introduced 
the vertex operators Z\^,(ff) defined in the regions 1^ according to Eq. (|4.2|l . We assume that these vertex 
operators can be analytically continued to the region (|i^'.l|l . It means that the spectrum in the regime II 
consists of three quadruplets of fundamental particles z*^, , i G Z3 (e and e' are eigenvalues of the operators 
Qi+i and Qi+2 respectively) with the masses 

„...M„.„f, M„ = iir(a)r(|)r(|). (f.4, 

and their bound states. The S matrices of the fundamental particles are derived from the commutation 
relations of the vertex operators Z*^, (G) . From the free field representation described in the present paper 
we can find that 

K,e'S^l)K,e'S^2) = ~ ^P.+i(^l-^2)???:5p.+.(01-e2):?:p:^e^ (02)^^34(^1), (F.5a) 

Z^,, (00^:^,(02) = ee" ^ S,Ue, - e,)%lZ^^l,,{B,)Zl,^{e,) (F.Sb) 
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with 

dt tanh ^ sin Qt 



- , - , , ( [°° dt tanh ^ sin 

s,«,)K = s,«,):: = „p(-.y^ 



2n+2 I 1 
7rp ' p 2 



T— r \ 7rp p ^ 1 V "^^^ V ^ I \ P ^ } 

-^-^ F2 ie 1 2k+1 I if-LSnilAr/ iS, 2ra+2 , 1 V (^■^) 

"=0 ^ p ^2)^\ 1^^~^2)^\ If^^ p ^2) 



cosh - „ „ . sin - 

Sp{Q)X~- = Sp{9)zX = - ^^^^ iJ'_g -Sp{9)XX, Sp{9)lt = Sp{9)t+ = - ^^^^ i^^g -Sp{9)XX- 

p p 

Note that the matrix Sp{9) can be expressed in terms of the sine-Gordon S matrix: 

5p(0) =itanh(f + ifi)5p(0 + if ) . ESI) 

The proof of the relations l|F.5b(l . I|F.6(I is similar to that of the commutation relations for Z^{9) given in 
Appendix ICl and we omit it. 

The interpretation of the result is the following. Consider the S matrix of two particles z* with the rapidity 
6*1 and with the rapidity 02- Let 9 ~ 9i — 92- The physical scattering process takes place for 9 > 0. Then 
for j = i the 5* matrix is given by the tensor product Sii{9) = —Sp-^-^{9) (8) Sp.^^i^)^ where the first tensor 
component acts on the space related to the topological charge Qi+i and the second one on the space related 
to the topological charge Qi+2- For j ~ i + 1 the S matrix is given by S'i.i+i(6') — QiQi^iSp.^^i^)^ while 
for j = i — 1 the S matrix is given by Si+i^i{9) — QiQi+iSp'.^^{^9), where Sp.^^{9) acts in the space of the 
topological charge Qi+2- 

Now consider the bound states. There are two types of bound states. Bound states of the first type 
are related to the poles of the S matrix for two fundamental particles of the same kind and are completely 
analogous to those in the regime I, described by Eqs. (|2.11|) . (|2.12|) . Explicitly, for any pi < 1 {i = 1,2,3) 
we have series of poles of the S matrices Sjj{9) with j = i ± 1 at the points 9 — iui^n with given by 
Eq. H2.11|l . The masses of the corresponding bound states 6^'*'", v — +,S, A, — are given by 

Mj;i,„ = 2mj sin— j i, n=l,2,..., np^ < 1. (F.7) 

It is important to note that 

M,.,,i = M,.,- 1 = 2Mo sin ^ sin i. (F.8) 
It can be checked that the quadruplets W'^'^ and are the same. 
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Bound states of the second type correspond to the poles of the S matrix Sij{9) for i ^ j. Let k ^ 
Then the function Sk{0) possesses poles at the points 

^ in - ^-^Pk - inpkn, n = 0,1,2,..., n < — - i. (F.9) 
z pk z 

The squares of the masses of the corresponding bound states read 

ml. = (sin^ ^ + sin^ ^ " ^ - ^ Bin ^ cos + npku) ) . (F.IO) 

We denote these bound states as z^^", s,e' = ±. Their multiplet structure is the same as that of the 
fundamental particle z''. Moreover, 

rrikfi = Mo sm = mk, 

and the lightest bound states z^^? coincide with the fundamental particles z^^,. The other bound state 
quadruplets z'^'", n > 1 appear if pk < 2/3. Evidently, the masses of these states are in the range 

1 1 

ruk < rnk,n < rrii + rrij, n = l,2,..., n< -. 

Pk ^ 

We suppose to give a detailed description of the field theory corresponding to the regime II in another 
publication. 
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